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The new method of nonperturbative calculation of the beta function in the 
lattice gauge theory is proposed. The method is based on the finite size scaling 
hypothesis. 
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! Monte Carlo simulations of SU(2) lattice gauge theory have shown that there are 

O I rather large deviations from asymptotic scaling behaviour in the range of coupling 

^ ' constants accessible with today's computing power. This raised the question of whether 

O ■ we are able to see continuum physics in these simulations. The answer requires the 



knowledge of the /5-function away from the asymptotic regime where it is dominated 



! by the two leading terms in its perturbative expansion (gf-coupling constant) 

(3f{g) = -bo9'-hg'-..., 
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where Bq = -, bi 
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/\ • During the last years a large effort has been put into the numerical determination of the 

■ /3-function for the lattice gauge theories by measuring the deviation from the two-loop 

/3-function (1) [1,2]. 

We propose a new method of determination of the nonperturbative /3-function, 
which is based on the finite size scaling hypothesis. 

We consider SU(2) gauge theory at finite temperature on A^^ x A^,- lattices with the 
standard Wilson action 

^(^) = 4E(l-^Trt/A (2) 

where Up is the product of link operators around the plaquette. The number of lattice 
points in the space (time) direction N„{Nt-) and the lattice spacing a fix the volume 
and temperature 

V = (iV,a)^ T = l/(iV,a). (3) 
The /3-function is defined by the expression 

«.)^-.^. (4, 
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Then for lattice spacing a one can obtain 



where is the renormahzation group invariant parameter. Since in Monte Carlo 
simulations the thermodynamic functions are calculated in the units of lattice spacing 
a, formula (5) determines the temperature dependence of these functions according (3). 
In the asymptotically free (AF) regime (1) we obtain the well known formula 



„A. = i;to^) = exp{-i|ln(M^)-^}, (6) 

which is valid in the region < 1. 

On the other hand, SU(2) gauge system at finite Nr and undergoes the deconfinement 
phase transition at g^^l. The new nonperturbative method for the calculation of 
/3-function is needed, which is not connected with the expansion (1). 

Our approach is based on the two points: i) translation into a more conventional 
statistical mechanical definition of /3-function and ii) the finite size scaling theory and 
phenomenological renormahzation. As in the standard spin systems, let us make the 
infinitesimal transformation of the lattice spacing a — > a' = 6a = (1 + Ab)a. Then 



- = - hmf a . ^M:i^U - hm ^ = - hm t . 
da 6^1 V ba- a ) Ab 6-»i db 

We obtain the new definition of /3-function for SU(2) lattice gauge system 



(7) 



It has been shown in the finite size scaling theory that on the finite lattice x 
{Nt- fixed) the order parameter (L), the susceptibility x correlation length ^ 

can be expressed in the following form (see, for example [3]) 

0{g-^N,)^NjQo{g-\K). (9) 

Here O represents (L), x and ^, cu — — /9,7, z/ is the corresponding critical index. 
Scaling function Q has some special dependence on g~^ and N^^, but this is out of our 
consideration. 
For example 

(L) = N;^Q^l){9-\N,), (10) 

The existence of the scaling function Q allows to develop a procedure to renormalize 
the coupling constant g^^ by the use of two different lattice sizes N^- and N^. Let us 
fix the physical size L = N^a and make a scale transformation 

{a ^ a' = ba, , . 
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Then the phenomenological renormahzation is defined by the equation 



Q{g-\N,)^Q{{g')-\Njh) 



(12) 



It expresses that the scahng function Q remains unchanged if the lattice size is rescaled 
by a factor h and the inverse couphng g~^ is shifted to {g')~^ simultaneously. Taking 
the derivative with respect to the scale parameter b of the both sides of (12) and using 
(8) it is easy to obtain the expression 



dg~^ d\nQ/d\nN^ 



da din Q/dg~'^ 



(13) 



The approximation of the derivative with respect to N^^ by the finite difference yields 
the final formula for the /3-function 



dg 
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hi 
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'N' 



da 



dQ{N,) dQiK) 
dg~'^ dg~'^ 



Q{N,)Q{K) 



1/2 
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So far we have considered the scaling function Q. It is very interesting to apply the 
same analysis to the correlation length ^. Using (10) we obtain instead of (12) 



i{g,K) ag',K) 



NL 
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In the case of the large enough lattice size [N^. — > oo) the dependence of the correlation 
length ^ on N„ become negligible and (15) yields 



i{g)^hi{g'). 



(16) 



This is the renormahzation group equation for the bulk system, which is wellknown in 
the standard spin theory. From this by the manipulations mentioned above one can 
obtain the expression for the /9-function 



dg 



(17) 



Substituting (17) into (5) we have an extremely simple formula for the lattice spacing 
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The first attempt to calculate the correlation length C,{g) has been made in [4], but with 
rather pure statistics and on rather small lattices {N„ = 18, = 3, 4, 5). We can see 
that the new method of calculation of the nonperturbative /5-function for SU(2) lattice 
gauge thermodynamics only needs the correct calculation of the correlation length ^ in 
the wide coupling constant interval. 

On the other hand, the best studied quantities in MC lattice calculations of SU{N) 
gauge theories are the string tension -^/a and the deconfinement transition temperature 
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Tc- The MC calculations give the values of the critical coupling j3^'~' of the deconfine- 
ment transition and dimensionless string tension {^/a)MC- The values (5^'-' were found 
for the finite lattices and the extrapolation to spatially infinite volume ("thermody- 
namical limit") A^^- — > oo has been done (see Ref. [5] and references therein). For 
the SU (2) gauge theory the MC values of the critical couplings 13^^ are presented in 
Table 1 for different Nr. One observes a rather strong dependence of Tc/Af^ on Nr. 
This means that the perturbative AF relation (6) does not work even on the largest 
available lattices. This fact is known as an absence of the asymptotic scaling. 

It has been proposed in Ref. [2] that a deviation from the asymptotic scaling can 
be described by a universal non-perturbative (NP) beta function, i.e. I3f{g) is the same 
for all lattice observables and it does not depend on the lattice size if A^^^, A^^- are not 
too small. The following ansatz was suggested [2]: 

aAr = A(/)i?(/), (19) 

where R{g^) is given by Eq. (6) and \{g'^) is thought to describe a deviation from 
the perturbative behaviour. The equation (6) has been expected at — so that 
an additional constraint, A(0) = 1, has been assumed. The values of Tc/A^^ can be 
calculated then as 

A simple formula for the function X{g'^) was suggested [2]: 

Parameter C3 in (21) and a new one, T*/A^^ = const, were considered as free pa- 
rameters and determined from fitting the MC values of T^/A^^ at different A^^- to the 
constant value T*/A^^. This procedure gives: 

r;/A^^ = 21.45(14), C3 = 5.529(63) • 10"^ (22) 

In comparison to Tg/Af^ the much weaker A^^ dependent values of T^/ Af^ have 
been obtained, which become now close to the constant value T*/A^^ (22). 

In spite of the phenomcnological success of the above procedure of [2] the crusial 
question regarding the existence of the universal NP beta function, with docs not 
depend on the lattice size, is not solved and remains just a postulate. To answer this 
question we reanalyze the same MC data using the different strategy as in Ref. [2]. 
A principal difference of our analysis is that we do not assume the existence of the 
universal beta function and take into account the finite size effects of the lattice. 

Usually finite size scaling (FSS) in the vicinity of a finite-temperature phase transi- 
tion is discussed for lattice SU(N) gauge models, without trying to make contact with 
the continuum limit, i.e. the scaling properties are studied on lattices of A^^ x A^^ with 
fixed Nr and varying A^,^, and the model is viewed as a 3-dimensional spin system. In 
the continuum limit the FSS properties of these non-abelian models should, of course, 
be discussed in terms of the physical volume V = and the temperature T in the 
vicinity of the deconfinement transition temperature Tc. We will study here how the 
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scaling behaviour of the continuum theory emerges from the lattice free energy on 
arbitrary lattices, i.e. when varying Nj- and A^^-- 

On a lattice N,,- x A^^ the length scale L and the temperature T are given in units 
of the lattice spacing a, therefore it is advantageous to replace the length scale L by 
the dimensionless combination 

TV 

LT = ^ (23) 

Using this ratio the singular part of the free energy density is described by a universal 
finite-size scaling function [3,5] 

where /9, 7, v are the critical indexes of the theory, the scaling function Q depends 
on the reduced temperature t = {T — Tc)/Tc and the external field strength h. 

Next we consider y — N^/Nt fixed, varying A^^- and therefore A^,- accordingly as 
is needed to reach the continuum limit. Rescaling Nfj and A^^ by a factor h leads to 
a phenomenological renormalization g'{g,b,y) by the following identity for a scaling 
function Q 

■(!)"") iV./^.(^)""'). (25) 

It follows from (25) 

t{g,N,)^t{g\Nr/b). (26) 

In general the reduced temperature t — {T — Tc)/Tc is a comphcated function of 
the coupling (3 = 2N/g'^, which in the visinity of the critical temperature Tc can be 
approximated by [5] 



1 - 



(27) 



This approximation reproduces the correct reduced temperature in the continuum limit, 
which is easy verified by using (4). Taking the derivative with respect to the scale 
parameter b of the both sides of (26) and using (8) and (27) it is easy to obtain the 
expression for the beta function: 

Pf{g) = -Bo{N,)g' - B,{K)g\ (28) 

where 

m\ boPcJdlnNr ^29) 

Then the equation (4) leads to 

a^L = exp(-^) (5o/)-^^/^^°. (30) 

Using (3) and (30) one can obtain the critical temperature Tc. The only problem is 
the calculation of the derivative dPc/dlnNr in expression (29). For the SU{2) gauge 
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theory the calculation has been made by fitting the MC data for the critical couplings 
P^*^ with a Sphne interpolation and numerical differentiation of this curve, The result 
of the calculation is represent in Table 1. In comparison to T^/A^^ the much weaker 
dependence on Nt of the critical temperature T^/Al is observed. 
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Table 1: MC data for Pc are taken from [5]. The values of TjAf^^ are calculated from 
(6). Our results for TJAl are obtained from (30). 
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TjAr 


dPc 

(INr 


TJAl 


2 


1.880 


29.7 






3 


2.177 


41.4 


0.158 


25.22 


4 


2.299 


42.1 


0.086 


25.46 


5 


2.373 


40.6 


0.063 


25.38 


6 


2.427 


38.7 


0.045 


24.13 


8 


2.512 


36.0 


0.040 


24.24 


16 


2.739 


32.0 


0.017 
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